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■ Abstract 

. Top quark polarization in e~^e~ annihilation into tt is calculated for linearly 

\^ I polarized beams. The Green function formalism is applied to this reaction near 

0^ ' threshold. The Lippmann-Schwinger equations for the S-wave and P-wave 

r"| ' Green functions are solved numerically for the QCD chromostatic potential 

given by the two-loop formula for large momentum transfer and Richardson's 
QhI ansatz for intermediate and small momenta. S-P-wave interference contributes 

^ ■ to all components of the top quark polarization vector. Rescattering of the 

decay products is considered. The mean values {ni) of the charged lepton four- 
momentum projections on appropriately chosen directions n in semileptonic 
^ ■ top decays are proposed as experimentally observable quantities sensitive to 

top quark polarization. The results for {n£) are obtained including S-P-wave 
interference and rescattering of the decay products. It is demonstrated that 
for the longitudinally polarized electron beam a highly polarized sample of top 
quarks can be produced. 
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1 Introduction 



The top quark is the heaviest known elementary object. Precision studies of its inter- 
actions will lead to profound progress in particle physics and may open an exciting 
new window towards the very high mass scale. Its large mass allows to probe deeply 
into the QCD potential which governs the dynamics of the nonrelativistic tt system. 
Such a system will provide a unique opportunity for a variety of novel QCD studies. 
It is likely that precise studies of top quark production and decays will reveal new 
important information on the mechanism of electroweak symmetry breaking. The 
analysis of polarized top quarks and their decays has recently attracted considerable 
attention |1|, |^. For nonrelativistic top quarks the polarization studies are free from 
hadronization ambiguities. This is due to the short lifetime of the top quark which 
competes favorably with the formation time of top mesons and toponium resonances. 
Therefore top decays interrupt the process of hadronization at an early stage and 
practically eliminate associated non-perturbative effects. 

Threshold production of top quarks at a future electron-positron collider will al- 
low to study their properties with extremely high precision. The dynamics of the top 
quark is strongly influenced by its large width ^ 1.5 GeV. Individual quarkonium 
resonances can no longer be resolved, hadronization effects are irrelevant and an ef- 
fective cutoff of the large distance (small momentum) part of the hadronic interaction 
is introduced 0, P]. This in turn allows to measure the short distance part of the 
potential, leading to a precise determination of the strong coupling constant [0 . The 
analysis of the total cross section combined with the top quark momentum distribu- 
tion will determine its mass rrit with an accuracy of at least 300 MeV and its width 
Tt to about 10%. For a Higgs boson mass of order 100 GeV even the tiH Yukawa 
coupling could be indirectly deduced from its contribution to the vertex correction 
[Q. Additional constraints on these parameters can be derived from the forward- 
backward asymmetry of top quarks and from measurements of the top quark spin. 
Close to threshold, for E = ^/s — 2mt <C m^, the total cross section and similarly the 
momentum distribution of the quarks are essentially governed by the 5'-wave ampli- 



tude, with P- waves suppressed ~ ~ yi?^ + F^/mj ^ 10~^. The forward-backward 
asymmetry and, likewise, the transverse component of the top quark spin originate 
from the interference between S- and P-wave amplitudes and are, therefore, of order 
(3 ~ 10~^ even close to threshold. Note that the expectation value of the momen- 
tum is always different from zero as a consequence of the large top width and the 
uncertainty principle, even for P = 0. 

It has been demonstrated 0] that the Green function technique is particularly suited 
to calculate the total cross section in the threshold region. The method has been 
extended to predict the top quark momentum distribution. Independent approaches 
have been developed for solving the Schrodinger equation in position space and the 
Lippmann-Schwinger equation in momentum space ||l^, |ll| . The results of these two 



methods agree very well. One of the most important future applications will be the 
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determination of rrit and c^s [§• A further generalization then leads to the inclusion of 
P-waves and, as a consequence, allows to predict the forward-backward asymmetry 



12]. It has been shown |13] that the same function (p^{p,E) which results from 
the S-P-wave interference governs the dynamical behaviour of the forward-backward 
asymmetry as well as the angular dependence of the transverse part of the top quark 
polarization. The close relation between this result and the tree level prediction, 
expanded up to linear terms in (3, has been emphasized. The relative importance of Z 
versus 7 and of axial versus vector couplings depend on the electron (and/or positron) 
beam polarization. All predictions can, therefore, be further tested by exploiting their 
dependence on beam polarization. In fact the reaction e~^e~ ti with longitudinally 
polarized beams is the most efficient and flexible source of polarized top quarks. At 
the same time the longitudinal polarization of the electron beam is an obvious option 
for a future linear collider. 

In the present article the polarization dependent angular and momentum distribu- 
tions of top quarks are studied (neglecting effects of CP violation, the corresponding 
distributions for i antiquark can be obtained through a CP transformation) and the 
results of |13| are expanded in three directions: 



1. Normal polarization. The calculation of the polarization normal to the produc- 
tion plane is a straightforward extension of the previous work [|I3|. It is based 
on the same nonrelativistic Green function as before, involving, however, the 
imaginary part of the interference term ip^{p, E). A component of the top quark 
polarization normal to the production plane may also be induced by time rever- 
sal odd components of the 7tt- or Ztt-coupling with an electric dipole moment 
as most prominent example. Such an effect would be a clear signal for physics 
beyond the standard model. The relative sign of particle versus antiparticle po- 
larizations is opposite for the QCD-induced and the T-odd terms respectively 
which allows to discriminate between the two effects. Nevertheless it is clear 
that a complete understanding of the QCD-induced component is mandatory 
for a convincing analysis of the T-odd contribution. 

2. Rescattering of decay products. Both t quark and i antiquark are unstable and 
decay into W~^b and W~b, respectively. Neither b nor b can be considered as 
freely propagating particles. Rescattering in the tb and bi systems affects not 
only the momenta of the decay products but also the polarization of the top 
quark. Moreover, in the latter case, when the top quark decays first and its 
colored decay product b is rescattered in a Coulomb-like chromostatic potential 
of the spectator t, the top polarization is not a well defined quantity. Instead 
one can consider other observables, like the angular momentum of the Wb sub- 
system, which coincides with the spin of top quark if rescattering is absent. 
Rescattering corrections are suppressed by ctg. The resulting modifications of 
the momentum distribution are, therefore, minor and as far as the total cross 
section is concerned can even be shown to vanish [11^ [THl. In contrast the 
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forward-backward asymmetry as well as the transverse and normal parts of the 
top quark spin are suppressed by a factor ~ (3. Hence they are relatively more 
sensitive towards rescattering corrections. The treatment of rescattering fol- 
lows to some extent the formalism of |P, |12|. However, instead of the Coulomb 
potential we employ the full QCD potential. 

The lifetime of the ti system is about (2rt)^^. At the time of the decay of 
the first constituent the typical potential energy of tb (or bt) is about rrit ~ 
agy/mtTt where the approximate identity between the two quantities is a nu- 
merical coincidence valid for rrit of about 180 GeV. At the time of the second 
(spectator) decay the mean separation of b and b has increased to Tf and the 
corresponding potential energy is reduced to as Tt. Thus, rescattering in the bb 
system is less important and will be neglected. 

3. Moments of the lepton angular distribution. The direction of the charged lep- 
ton in semileptonic decays is optimally suited to analyze the polarization of the 
top quark 0. The reason is |jl8[ that in the top quark rest frame the double 



differential energy-angular distribution of the charged lepton is a product of the 
energy and the angular dependent factors. The angular dependence is of the 
form (1 -|- Pcos6), where P denotes the top quark polarization and 6 is the 
angle between the polarization three-vector and the direction of the charged 
lepton. Gluon radiation and virtual corrections in the top quark decay practi- 



cally do not affect these welcome properties |]19|. It is therefore quite natural 
to perform polarization studies by measuring the inclusive distributions of say 
/x+ in the process e+e" — >■ t{^~^i>^b)t{jets). This can be also convenient from 
the experimental point of view because there is no missing energy-momentum 
for the t subsystem. From the theoretical point of view the direction of the 
charged lepton is advantageous since it is equivalent to the top quark polariza- 
tion when rescattering is absent and remains defined even after bi rescattering 
is included. However, the semi-analytic calculation of this effect is difficult be- 
cause production and decay mechanisms are coupled. In this article moments of 
experimentally measurable distributions are proposed which are less difficult to 
calculate but still retain the important information on top quark polarization. 
These moments, the mean values (ni) of the charged lepton four-momentum 
projections in semileptonic top decays are calculated including S'-P-wave inter- 
ference and 0{as) rescattering corrections. 

The outline of this paper is as follows: In Sect. |] the relations between the mo- 
mentum, angular and spin distributions of top quarks and S- and P-wave Green 
functions are established. The formalism is developed in sufficient generality such 
that the application to other reactions like 77-fusion into ti is obvious. Simple rules 
for the translation from tree level approximation to the case including the QCD po- 
tential are presented. This section expands and details our previous results, as far as 
the polarization in the production plane is concerned, and presents new results for 
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the normal polarization. In Sect. ^ a formalism for the rescattering of the top quark 
decay products is presented using the spin projection operators for the top quark. 
The influence of the rescattering on the various differential (momentum, angular, po- 
larization) distributions is studied and limitations of the formalism in the case of bt 
rescattering are discussed. In Sect. |^ the moments {n£) are defined and calculated. 
In Sect. 1^ the numerical solutions of the Lippmann-Schwinger integral equations are 
presented. Predictions are given for the longitudinal, perpendicular and normal po- 
larization including predictions for the corresponding moments (ni). The dependence 
on the various input parameters and the effects of rescattering are studied. Various 
technical aspects are relegated into the appendices. 

2 Green functions, angular distributions and quark 
polarization 

2.1 The nonrelativistic limit 

Top quark production in the threshold region is conveniently described by the Green 
function method which allows to introduce in a natural way the effects of the large 
top decay rate and avoids the summation of many overlapping resonances. The 
total cross section can be obtained from the imaginary part of the Green function 
G(x = 0, x' = 0, E) via the optical theorem. To predict the differential momen- 
tum distribution, however, the complete x-dependence of G(x, x' = 0,-E') (or, more 
precisely, its Fourier transform) is required. In a calculation with non-interacting 
quarks close to threshold the forward-backward asymmetry, the angular dependent 
term (~ cos^) of the longitudinal part and the transverse part of the top quark 
polarization are all proportional to the quark velocity (3 and originate from the in- 
terference of an S'-wave with a P-wave amplitude. These distributions are described 
by V ■ G(x, x', E) |x'=o or, equivalently, by the component of the Green function with 
angular momentum one. The connection between the relativistic treatment and the 
nonrelativistic Lippmann-Schwinger equation has been discussed in the literature 
|12[. The subsequent discussion follows these lines. It includes, however, also the 
spin degrees of freedom and is, furthermore, formulated sufficiently general such that 
it is immediately applicable to other reactions of interest. The main ingredient in 
the derivation of the nonrelativistic limit is the ladder approximation for the vertex 
function F^. This vertex function is the solution of the following integral equation 
(Fig. D: 

Tc = C + / ^ [-f^c^s) D,,{p - k)YSAk + l)Tc{k, q)S^{k - |)7^ (1) 

with C = 7^ or 7/^75 in the cases of interest for e+e^-annihilation. The conventions 
for the flow of momenta are illustrated in Fig. The four-momenta are related to 
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t 



t 

Figure 1: Lippmann-Schwinger equation ([^) in diagrammatical form. 

q/2 + p 



P 

q/2-p Figure 2: Definition of the four-momenta. 



the nonrelativistic variables by 

q = {2mt + E,0) 

P = (0,p) (2) 
k = {ko,k). 

In perturbation theory the ladder approximation is motivated by the observation that 
for each additional rung the energy denominator after loop integration compensates 
the coupling constant attached to the gluon propagator. This is demonstrated most 
easily in Coulomb gauge. Contributions from transverse gluons as well as those from 
other diagrams are suppressed by higher powers of /3 ~ as- The gluon propagator is 
thus replaced by the instantaneous nonrelativistic potential: 

4 

- - 4:7iasD^^{p) i V{p) 5^o5uo. (3) 

The dominant contribution to the integral originates from the region where |k| rrit. 
Including terms linear in k, quark and antiquark propagators are approximated by 
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SF{k + ^ 



— + /Co ^s::^ + *~ 



2 ' ^ 2mt ' 2 



^ 2" ~ '^O - 2^ + 



The "elementary" vertex C is independent of ko- (Within the present approximations 
this is even true if C does depend on k as is the case in the analogous treatment 
of 77 ti discussed below.) Up to and including order [3 terms a self-consistent 
solution of the integral equation ([^) can be obtained if Vc is taken independent of /cq 
and the nonrelativistic spins of t and i. The ko integration is then easily performed 
and the integral equation simplified to 

In the calculation of the cross section for the production of t plus i with momenta 
g/2 ±p and spins s± respectively traces Ti of the following structure will arise: 

H = Tr{P+(|+p,s+)rcP4|-P,S-)fc'}, (6) 

^ / x ^±mt 1 + 75^ . . 

= ^' 

where we allowed for mixed terms with C different from C arising e.g. from vector- 
axialvector interference. Expanding again up to terms linear in k, this trace can be 
transformed into 

n = TiiS+TcSJc'}, (8) 
S± = (9) 

S = 7757°= (o °), (10) 
with the nonrelativistic reduction defined through 

fc(P,£) = A,(i-E-2)r,(p,£)(i-E-I)A_. (11) 

It is thus sufficient to calculate the "reduced" vertex function Tq- Dropping again 
terms of order k^, the corresponding integral equation is cast into a particularly simple 
form _ 
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where C(p) is defined in analogy to rc{p,E) in (]Tl|). Consistent with the nonrela- 
tivistic approximation only the constant and the linear term in the Taylor expansion 
of the elementary vertex will be considered]^ 



C(p)=C(0) + D-p. 



(13) 



The matrices C(0) and D may in general depend on external momenta, polarization 
vectors or Lorentz indices. A self-consistent solution for the vertex Tc is then given 
by 

f c(p, E) = C(0)/Cs(p, E) + B-p /Cp(p, E). (14) 
The scalar vertex functions /Cs,p depend on the modulus of the three momentum 



P = IPI 

(which should be distinguished from the four-momentum p) and the nonrelativistic 
energy E only. They are solutions of the nonrelativistic integral equations 



ICs{v,E) 
/Cp(p,E) 



1 + 
1 + 



(^^p-k)^^^^ 
d^k p ■ k 



E-^ + iTt 

rrtt '' 



(27r)3 p2 



V{p - k) 



/Cp(k,E) 
E-^ + iVt 



(15) 
(16) 



and are closely related to the Green function ^(p,x, E') which, in turn, is a solution 
of the Lippmann-Schwinger equation 



^ - — + ^r^ 



6;(p,x,E) = e'P'' + 



(2^ 



y(p-k)e;(k,x,E). 



;i7) 



Let us denote the first two terms of the Taylor series with respect to x by G and F 
respectively: 

g{p, x, E) = ^(p, E) + X ■ zp F(p, E) + ... 
They are solutions of the integral equations 



with 



Go(p,i?) 



G(p,E) = Go(p,i?) + Go(p,i?)/^V(p-k)G(k,E) 
F(p,E) = Go(p, E) + Go(p, E) I ^^^^3 ^^^ V{p - k)F(k, E) 

1 



E-S1 + ^T^ 



In the notation of Ref.[|§ one gets: C(0) = A+OqA- and D = A+ - 2^{Oo, 7}+ + O 



(18) 
(19) 

(20) 
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and the relation between Green function and vertex function 

^(p, E) = Go(p, E)ICsiv, E), F(p, E) = Go(p, E)}Cp{p, E) (21) 

is evident. In the case of e"'"e~-annihilation top production proceeds through the 
space components of the vector and axial vector current. The relevant elementary 
vertex C(p) is given by 

T^(P) = A+7,A- (22) 
77y5(p) = A+(— )(7xp),A_ (23) 

TTlt 

for vector and axial current respectively. Production of tt in 77-fusion would lead to 
an elementary vertex of the form 

C(0) oc i{ei X €2) ■ ne-A+75A- (24) 
1 

D oc — A+ [(ei ■ e2)(ne- ■ 7)ne- + (e2 ■ 7) ei + (ei ■ 7) €2] A_, (25) 

nit 

with ei, £2 the polarization vectors of the photons, and the present formalism applies 
equally well. This case has been studied in |^ . 



2.2 Cross sections and distributions 

In the nonrelativistic limit (and defining Et = — rrit) the two particle phase space 
for the production of a stable t- and t-quark is conveniently written in the form 



1 n r-. 1 



-o27r5(p° - + p?) -^2n6{pl - ^mj + p?) = (26) 

APl Zp^ 

= ^ /^^"/^2vr5(E,- J^)2vr5(i?,-^)2vr5(E-E,-E,). 
4mfy 271 27r (27r)3 ^ 2m/ ^ * 2m/ ^ *^ 

The first two delta functions represent the on-shell conditions for top and antitop 
quarks. In the case of unstable particles they are replaced by Breit-Wigner functions 

2vr5(E - ;^) -— — — , (27) 



2^* iE-& + {rt/2y 



see also discussion in Sect. |^.l . The dEt integration is trivial. The dEt integral can 
be performed explicitly with the help of the residue theorem. For this second step it 
is essential that the matrix element is independent of Et for fixed E, an assumption 
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evidently fulfilled by the amplitude derived above. The phase space integration is 
thus cast into the following form 

Therefore only the combinations 

iGo/Csl' = IGI^ (29) 
/CsGo/C^G* = GF* (30) 

appear in the calculation of the ^-wave dominated differential cross sections and the 
5'-P-wave interference terms. 



2.3 Top production in electron positron annihilation 



With these ingredients it is straightforward to calculate the differential momentum 
distribution and the polarization of top quarks produced in electron positron annihi- 
lation. Let us introduce the following conventions for the fermion couplings 



Vf = 21 j — Aqf sin^ 9t^, a/ 
P± denotes the longitudinal electron/positron polarization and 



2ll 



1 - P. P_ 



(31) 



(32) 



can be interpreted as effective longitudinal polarization of the virtual intermediate 
photon or Z boson. The following abbreviations will be useful below: 



03 

04 
d 



= ?e?t + (^e + «e)^t d + 2q^qtVeVtd 

= 2veaeVfd'^ + 2q^qtaeVtd 

= AvetteVtatd^ + 2qeqtaeatd 

= 2{vl + al)vtatd'^ + 2qeqtVeatd 
1 s 



(33) 



16 sin 6'w cos^ 6'w s 



Ml 



The differential cross section, summed over polarizations of quarks and including 
5"- wave and S'-P-interference contributions, is thus given by 



dp^ 



1 - P.P. 



+(03 + X^i) ( 1 



12a. 



X02) ( 1 - 



P 



IGcts 



\G{V,E)\' + 



StT ) TTlt 



Re{G{p,E)F*{p, E)) cosi? 



(34) 



The vertex corrections from hard gluon exchange for S'-wave [53] and P-wave |2| 



amphtudes are included in this formula. It leads to the following forward-backward 
asymmetry 

AMv,E) = CFBix)MP^E), (35) 



with 



(p^ = Re ip, and 



1 Q3 + XQ4 .ORN 

<^fb(x) = ■ , (36) 

2 ai + xci'2 



(l-4^s/37r) p FXp,E) 



This result is still differential in the top quark momentum. Replacing E) by 



(i-wjr<ipgF-(p.g)G(p.^) 

' (l-8a,/3,r) /„»" <ipp2 |G(p, ' 

one obtains the integrated forward-backward asymmetry]^. The cutoff p^, must be 
introduced to eliminate the logarithmic divergence of the integral. For free particles 
(or sufficiently far above threshold) one finds for example 

$(E) = xl^+ ^ In ^ . (39) 
V mt y/mtEn nit 

This logarithmic divergence is a consequence of the fact that the nonrelativistic ap- 
proximation is used outside its range of validity. One may either choose a cutoff of or- 



der rrit or replace the nonrelativistic phase space element pdp/rrit hj pdp/ymf + p^. 
In practical applications a cutoff will be introduced by the experimental procedure 
used to define tt-events. 



2.4 Polarization 

3. t 



Figure 3: Definition of the spin 
directions. The normal compo- 
nent sn points out of the plane. 



•^For the case without beam polarization this coincides with the earUer result [Q, as far as 
the Green function is concerned. It differs, however, in the correction originating from hard gluon 
exchange. 
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To describe top quark production in the threshold region it is convenient to ahgn 
the reference system with the beam direction (Fig. ^) and to define 



Sn 



X nt 



Up 



X 



(40) 



Sn X sii 



In the hmit of small (3 the quark spin is essentially aligned with the beam direction 
apart from small corrections proportional to (3, which depend on the production angle. 



A system of reference with sy defined with respect to the top quark momentum |25 
is convenient in the high energy limit but evidently becomes less convenient close to 
threshold. 

The differential cross section for production of a top quark of three-momentum p and 
spin projection +1/2 on direction s+ rea^ 



d^(^{p, Sh 



fa 



(41) 



The three-vector 7-* characterizes the polarization of the top quark. Its components 
can be obtained from the formula 



(i^a(p,s+) d^a{p,s 



dp"^ 



dp^ 



d'a{p,sf) rfV(p,s^O 



dp^ 



+ 



dp^ 



(42) 



where s 



s , and sn- Including the QCD potential one obtains for the three 



components of the polarization 
^||(P,^,X) 



qi(x) + ql(x)v^K(p,^) cos^ 

C±(x)</?j,(p,^) sin 7? 
Cn(x)v5i(p,^) sin ■(9, 



(43) 
(44) 
(45) 



c^ix) = clix) = (1 - X') 

C^{X) = -7^ —— , Cn X = -- —— = -Cfb X , 



with = ImLf, and ip{p, E) as defined in (|37|) . The momentum integrated quantities 
are obtained by the replacement <y9(p, E) $(i?). The case of non-interacting stable 
quarks is recovered by the replacement an obvious consequence of (^). 

Let us emphasize the main qualitative features of the result: 

*In the non-relativistic approximation S-|_ is the same in the top quark rest frame and in the ti 
center-of-mass frame up to corrections of order 0^ which are neglected. 
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Top quarks in the threshold region are highly polarized. Even for unpolarized 
beams the longitudinal polarization amounts to about —0.41 and reaches ±1 for 
fully polarized electron beams. This later feature is of purely kinematical origin 
and independent of the structure of top quark couplings. Precision studies of 
polarized top decays are therefore feasible. 

Corrections to this idealized picture arise from the small admixture of P-waves. 
The transverse and the normal components of the polarization are of order 10%. 
The angular dependent part of the parallel polarization is even more suppressed. 
Moreover, as a consequence of the angular dependence its contribution vanishes 
upon angular integration. 

The QCD dynamics is solely contained in the functions or $ which is the 
same for the angular distribution and the various components of the polariza- 
tion. (However, this "universality" is affected by the rescattering corrections as 
discussed in Sect. |[) These functions which evidently depend on QCD dynam- 
ics can thus be studied in a variety of ways. 



• The relative importance of P-waves increases with energy, $ ~ ^Ejmt. This 
is expected from the close analogy between <I>r = Re $ and /?. In fact, the order 
of magnitude of the various components of the polarization above, but close to 
threshold, can be estimated by replacing $r — > p/m^. 

The Ci are displayed in Fig. ^ as functions of the variable x- Fo^' the weak mixing 
angle a value sin^6'w = 0.2317 is adopted, for the top mass mt = 180 GeV. As 
discussed before, assumes its maximal value ±1 for x = -pl and the coefficient C|| 
is small throughout. The coefficient C_l varies between +0.7 and —0.5 whereas Cn is 
typically around —0.5. The dynamical factors $ are around 0.1 or larger, such that 
the P-wave induced effects should be observable experimentally. 
The normal component of the polarization which is proportional to has been pre- 
dicted for stable quarks in the framework of perturbative QCD |2^, |2^. In the thresh- 
old region the phase can be traced to the tt rescattering by the QCD potential. For 
stable quarks, assuming a pure Coulomb potential V = —4as/3r, the nonrelativistic 
problem can be solved analytically and one finds 



hm [E-^^+zr)jG{p,E) = exp(^)r{l + zk/p) (47) 



lirn, [E-^ + tTt]F{p,E) = (l - t^) exp T{1 + ik/p), (48) 



with k = 2mtas/3 and hence 



2 1 - 4as/37r , , 
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Figure 4: The coefficients (|6]) for = = 180 GeV. 



The component of the polarization normal to the production plane is thus approxi- 
mately independent of E and essentially measures the strong coupling constant. In 
fact one can argue that this is a unique way to get a handle on the scattering of heavy 
quarks through the QCD potential. 

3 Rescattering corrections to the differential cross 
section 

3.1 Introductory remarks 

Close to threshold, the velocity /? of the top quark is of the same order of magnitude 
as the strong coupling. To include only the contributions to the differential cross 
section and polarization and discard rescattering corrections is thus inconsistent. To 
obtain reliable predictions for the observables it is mandatory to calculate the 0{a^ 
corrections. Some of these can be implemented without much effort as they are 
known since several years, and in fact they have already been taken into account in 
the preceding formulae: 

• the modification of the ttV and ttA vertices due to transverse gluon exchange 



for small /5 simply reduces to a multiplication by constant form factors |^0[ : 



G(p,E) — . (i-^)g(p,E) (51) 



13 



F{p,E) 



l-^)F(p,E). 



(52) 



This explains the corresponding factors in ([3^ ) and in the definition of V9(p, E) 
and ^{E) (|7D. 



Corrections to t or t decay can be taken into account by using the full 1-loop 
result for Ft 



If the total cross section were the only quantity of interest, all next-to-leading 
contributions would be included this way, because it is well known since quite some 
time that the excitation curve is free of (other) 0{as) corrections |31|. This result has 
also been re-established recently |T^, 0. But this statement does not hold for the 
differential cross section or polarization, since these are affected through the diagrams 
shown in Figs, ^a and Pd (and a third diagram with gluon exchange between the two 
6-quarks, which can, however, be neglected), i.e. corrections caused by interactions 
between the decay products of one of the top quarks with the second top. They will 
also be called rescattering corrections. 





a) 



b) 



Figure 5: Lowest order rescattering diagrams. 

Their calculation obviously requires a modification of the approach as at least the 
decay of one of the two top (anti) quarks has to be included. In this section a formalism 
will be pursued which attempts to retain the notion of quark polarization in the 
presence of rescattering of decay products. In Sect. ^ a slightly different strategy will 
be adopted which aims directly at predictions for experimentally observable lepton 
distributions, viz. appropriately chosen moments. Let us start with the first approach 
and write the differential cross section in the following form^ : 

+0{p) terms. (53) 

The three-vector s_|_ describes the spin of top quark, which in this subsection is con- 
sidered as a well defined quantity. The leptonic tensors L*-' are defined through 

= SmlS'^ , = -Sm]ie'^ , i,j = 1..2 . (54) 



^with the 3-axis identical to the electron beam axis 
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As already discussed in the preceding section the terms of order j3 arise from S'-P-wave 
interference, i.e. from interference between the amphtudes corresponding to vector 
and axialvector couphngs ^7 and tiZ^. In the approximation of the present article 
contributions suppressed by as(3 are neglected. Hence it suffices to consider order 
as corrections to contributions involving only vector couplings and the corresponding 
vertex functions /Cs(p,-E). represents the leading (vector- vector) part of the 

hadronic tensor 

= JdPS2{t;bW^) JdPS2{i-bW~) E ^^■^] + c-c. (55) 

where Aii denotes the amplitude for V ^ ti ^ bbW~^W~ . The momenta of top and 
antitop are given hj t = b + = q/2 + p, i =b + W~ = q/2 — p and e± denote the 
polarization vectors of W^. 

This approach can be adopted in a straightforward way as long as only factorizable 
diagrams are considered where the top quark decay products do not rescatter. Di- 
agram Fig. however, as a truly non-factorizable contribution leads to difficulties 
when discussing the top quark polarization as will become clear below. Neverthe- 
less this formalism will be presented for pedagogical reasons before a more powerful 
approach will be introduced in Sect. |[ 

To lowest order in the nonrelativistic approximation the amplitude Aii is given 

by 

Mfis^) = -{^)\ib)Ml-l')^-^f^A^l^A-/M'l>ib) 

■G{p,E){Gt{p,E) + Gt{p,E)) (56) 
with the nonrelativistic propagators introduced in (4) (see also App. 

Sf{p ± |) = ^(a± - ^) • G,/f(p, E), 

and where 

-SF{p + q/2)rySF{p-q/2) ^ A+-f'A^}Cs{p, E)Gt{p, E)Gt{p, E) 

= A+YA-GiP,E)iGtip,E) + G-tip,E)) (57) 

has been used. The amplitude Aif'\sj^) depends on s+ which for a stable top quark 
would be interpreted as its polarization four-vector. For unstable top quark this in- 
terpretation is still valid if the decay products do not interact. In such a case the 
total angular momentum of the W~^b subsystem is conserved. Thus, in principle at 
least, the angular momentum state l^jm) of W^b can be inferred from the measure- 
ment of the spin-dependent angular distributions of W and b. The quantum number 
m denotes an eigenvalue of the angular momentum projection on direction s+. It is 
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evident that one can identify m with the projection of top quark spin if there is no 
interaction changing momenta and polarizations of W and b after the decay. 

Inserting the expression for into (^) and (0), summing over polarizations e+ 
and e_, and performing the {bW~^) and {bW~) phase space integrations with the help 
of 

^ / dFS,{t; bW^) Y^i^yi^(\- = - 7') (58) 

and the corresponding identity for the {\)W~^ system, one re-derives the results given 
in the preceding section. The inclusion of the axialvector couplings and calculation 
of the S'-P-wave interference effects can be done in an analogous way. 



3.2 Rescattering in the ib system 

The rescattering in the ih system is given by diagram |^ and leads to an additional 
amplitude of the form 

-M- = (^)7(0I^W ^+(1-7'M+7'^A+7.A-/-(1-7')5f(A:-P-&) 

^-v{b)-Gt{v.E)G%E)[Gi{k,E)^Gj{k,E))D^,{k-v)gl^ • (59) 

In this amplitude all the top quark propagators have been substituted by their non- 
relativistic limits 

S^(t)^%k^Gt{v.E). 

As a consequence the top quark couples to gluons via the projection A+7^A+, 
which implies that only instantenous Coulomb-like gluon exchange contributes to 
leading order. In the calculation of the 0{ol^ corrections to the differential cross 
sections 

- JdFS2{t;bW+) JdFS2{t;bW-) J2 {Mf\s^)M'f^^ + M^^ Mf\s^)) 

only the interference of Ai"^ with the lowest order amplitude Aif\s+) contributes. 
A more detailed description of the calculation of this contribution is presented in 
App. At this point it should suffice to mention that it essentially involves choosing 
the appropriate integration contours and using the residue theorem. The result can 
be written in the following compact form: 

where da^^^dp denotes the top quark momentum distribution in the lowest order 
approximation, 

P|l = C^^Mp,E) + K^^{p,E) cos ^ , (61) 
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Pl = Pn = 0, (62) 

K = and y = — ^. (6d) 

l + 2y mi 

The two functions ipi and ip-^^ are defined as the following convolutions: 



r d^k Gfk E) arctan 

M,,E) ^ 2.,„/p-^n|k-p|)J-i^^— I- (64) 

^^'^ T./i, ,xP-(k-p)G'(k,E) 



^2(p,i5) = 2/7^n|k-p|) 



(27r)3 ^' "'^|p||k-pPG(p,E) 

r I k I 

xfl — -arctan ) (65) 

V |k-p| ^ ^ ' 

V^Jp,i?) = Re^^{j>,E). (66) 

It is noteworthy that rescattering in the th system does not affect the transverse 
and normal components of the top quark polarization. The top quark width acts as 
an infrared cutoff: in the limit — » the real part of ip2 becomes infrared divergent. 

3.3 Rescattering in the {Wh)i system 

The corrections due to rescattering of the top quark decay products on the spectator i 
can be calculated in an analogous way. They result from the interference of diagram 
||b and the lowest order amplitude M.f\s+). However, for diagram it is not 
possible to identify the total angular momentum j of the W^h system and the spin 
of the top quark. In fact the latter quantity is not well defined and j can be different 
from 1/2. However, only for j = 1/2 and m = ±1/2 there are non-zero contributions 
to order as corrections 

Summing over polarizations and calculating the phase space integrals one obtains 
with 

P|l = Cf^i(p,E) + kVr(p,^)cos^9 (68) 
P± = Kip^{p,E)sm^ (69) 
Pn = C?\Kij^{p,E)sm^ (70) 



where the additional function 



^,{p,E) = lmi;2ip,E) (71) 
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appears. 

At first glance this result may be surprising. One might have thought that any 
interaction after top decay cannot influence top polarization. There is however a 
problem with such a simple picture, because the spin of an intermediate (virtual) 
particle is not a well defined quantity and thus no observable. What has been calcu- 
lated by using the spin projection operator should be interpreted as the (bW) total 
angular momentum, which coincides with top spin on Born level only. Clearly the 
angular momentum of Wb may well be influenced by final state interactions. From 
a practical point of view it will be quite difficult if not impossible to measure this 
quantity. 

It definitely makes more sense to consider the angular distribution of the charged 
leptons from the decay t W^b l^vfi instead. Rescattering in the {Wb)i system 
produces a contribution which is non-f actor izable, i.e. it is impossible to write the 
final expression in the form 

dV(e+e- ^ l+vibW-b) ^ d^a{e+e- ^ tt) 1 dVt^uuM 

dpMEidni dp^ Tt dEidVli ^ ^ 

which is possible for all the other diagrams. 



3.4 Momentum distribution and polarizations 

Adding the contributions discussed in the preceding subsections one can calculate 
the order as correction to the differential cross section summed over spins. The 
momentum distribution is affected by tpi only: 

da'"! (iff"") , , , 

It follows that at order the total cross section for tt production is unaffected by 
rescattering corrections [14, 1^. Indeed 

J-0^,\G{p,E)\'Mp,E) = O (74) 



because for any real function /(|k — p|) the integral 



d^p f d^k 

'=]] 
is real. 

The function '?/'2(p,-E) gives the order corrections to the forward-backward 
asymmetry and the polarization vector 'P of the bW system. It is interesting to 
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note similarities between these formulae and the expressions: 



rose . 

rcsc , 



^Kip^{p,E) sint? 
^/t-?/'j(p, E) sin '(9. 



K ip^ (p, E) cos-i? 



(75) 
(76) 
(77) 

(78) 



Let us close this section with two remarks. The results for the unpolarized case 

Tt 



have been first derived in and are in accordance with ours. Following IT 
one can consider also the diagram with gluon exchange between b and b. However, 
it has already been argued in the Introduction that this contribution is suppressed. 
Our conclusion is supported by the results of |jl5|, |16| for unpolarized beams. 



4 Moments of the lepton spectrum 

It has already been stated that the optimal way to experimentally determine top 
quark polarization is based on the analysis of the charged lepton distribution in the 
process 

e+e" ^tt^ blubW. 

A particularly convenient distribution is the charged lepton angular distribution, 
which for free top quark decay in its rest frame is given by the formula 

- = -(l + |P|cos») (79) 

where 6 denotes the angle between the lepton momentum and top polarization T'. It 
is quite straightforward to generalize this formula to the case of factorizable contri- 
butions to ti production, c.f. ([7^). In fact this problem can be solved by using the 
Lorentz transformation from the top quark rest frame to the center-of-mass frame 
for the ti system. However the diagram of Fig.^D leads to a non-factorizable con- 
tribution with production and decay mechanisms coupled in an intricate way. The 
corresponding expression for the angular distribution of the charged lepton involves a 
complicated integral over the phase space and the three-momentum transfer between 
b and t. In particular integration over the phase space of the top quark is difficult 
because the direction of the charged lepton in t — >■ bW — >■ blu is fixed and one can- 
not use Lorentz symmetry to perform integrations. We were unable to simplify this 
expression. It is more convenient to consider moments of the charged lepton distri- 
bution which are given by manifestly Lorentz invariant integrals over the three-body 
Lorentz invariant phase space PS3{t;bli>). Moreover, these moments are measurable 
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and contain important information on polarization. From the calculational point of 
view Lorentz symmetry simplifies the problem considerably and closed expressions 
can be derived even for non-factorizable contributions. 

Let / be the four-momentum of the charged lepton and n denote one of the fol- 
lowing four- vectors, which in the ti rest frame read, c.f. (|iO|): 

nil = (0,S||) (80) 

= (0,sx) (81) 

nN = (0,Sn) (82) 

n(o) = (1,0,0,0) (83) 

The moments are defined as the average values of the scalar products {nl). One can 
consider the moments (jil) for a fixed three-momentum of the top quark, 

or, integrating over the direction of the top quark, the moments 

jiE,da4n, \^ia,ia4E, f*^' 

or, finally the moments {{{nl))) integrated also over p. The calculation of these 
quantities is straightforward because the integral over the three-body phase space 
PSslt^blu) can be calculated analytically. The hadronic tensor for the differential 
cross section is given by 

H^^=J^J^J dPS^it-blu) J dPS2{t;bW-) J:MM] (86) 

with the appropriate amplitudes including semileptonic decays. The correspond- 
ing tensor needed for the determination of the moments will be called Hij-n and simply 
contains the additional factor (nl) under the integrals 

H^r,n = 1^1^ I dPS,it;blu) J dFS2{t;bW-) J2M^Mlinl) (87) 
Using the following decomposition of the three-body phase space 

dPSsit; blu) = dW^ I dPS2{t;bW) J dPS2{W;lu) 



one can reduce the calculation to a sequence of Lorentz invariant integrations over 
two-body phase spaces. A more detailed description of the calculation can be found 
in App. p. 
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Let us, for illustrative purpose, discuss the Born- and S-P-wave interference part. 
In this case one arrives at 



(n/) = BR(t ^ blu] 



1 + 2t/ + 3i/^ 
4(1 + 2y) 



(89) 



with y (and k) as defined in (|63D. In the tt center-of-mass frame the space components 
of V are given by (H3])-(H3) and 



nit 



(90) 
(91) 



This clearly demonstrates that for the factorizable contributions the moments con- 
tain information equivalent to those obtained by using the approach of Sect. 2 and 
Sect. 3. Integrating over the emission angle of the top quark one derives the following 
expression for the angular independent part of the longitudinal polarization 



(0) 



Cf = -3((n||0)/(Ko)/)) 



(92) 



Similarly the coefficient characterising normal polarization can be obtained from 

(cf.(ii)) 

C^lmip{p,E) = - — (KO) /((^(o)0) ■ (93) 



TC 



Let us now return to the case with rescattering included. It requires the introduc- 
tion of a third convolution function that only enters through diagram ^d: 



^3(P,^) 



Im 



pi 



p- (k-p) ^2 

|p||k- p| 



1 



(94) 



(27r)3 ^' "'^G(p,E) 
The complete result for the moments, correct up to (and including) C(/3, as) reads 

l + 2y + 3y^ 



{nl) 



BR{t ^ blu 
itn){l + — 



4(1 + 2y) 



+ 2y){l + 2y + 3y^) 



rrit ^ rrit 

H — -nV H -n5V 

3 3 



(95) 



where the complete differential cross section as given by 

da da'-^^ , , _x 1 + 2^fb(p, ^) cos?? 



dpdflp 
^fb(p, E) 



dp 



At: 



with 



Cfb ¥'r(p, E) + -^K ^/'j,(p, E) 
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(96) 
(97) 



should be used as normalization. The change in "polarization" caused by rescattering 
looks as follows: 



2 + 3y-5y^- 12y^ 
;i + 2y)(l + 2y + 3y2 



K 



Q\2 



(p, E) cost? 



5Vi 



4(1 + 2y + 3y2 
3(1 - 3y'^ 



(98) 



2(l + 2|/ + 3|/2)^« 
3(l-4y + 3y2 
8(l + 2y + 3y 

0, 6Vo = 0. 



t/'j^(p, ii^) sin?? 
Cf^3(p, ^) sin2t? 



(99) 
(100) 



An important consequence of these expressions - which obviously differ from those 
given in Sect. |5 - is that the normal component of the polarization vector is unaffected 
by rescattering corrections, and that ( P^ and (PB| ) remain valid. The additional term 
multiplying tn in ( pSf ) and the two tps dependent terms on the other hand seem to spoil 
the beauty of and the similarity with the Born result. But it had to be expected that 
the rescattering corrections would destroy such a simple form. Numerically however, 
because of the large mass of the top quark, it will turn out below that these three 
terms can be neglected. 



5 Numerical results 

The QCD potential used in the numerical analysis has been described in |]ll|. It 
depends on the strength of the couphng constant. If not stated otherwise aj^(M|) = 
0.125 is adopted. For the top mass a value of 180 GeV is used and, correspondingly, 
Tt =1.55 GeV [1^. The effects of a momentum dependent width on cross section 
and momentum distributions have been studied in |11| and shown to be well under 
control and small. 

Step by step the predictions will be presented for the momentum distribution, the 
forward-backward asymmetry and the top quark polarization and thus for observables 
of increasing degree of complication. They will be shown to involve increasingly 
complicated rescattering corrections. For each of them the prediction based on S-P- 
wave interference will be given which allows to understand the qualitative features of 
the result. Subsequently the modifications due to rescattering will be introduced. 



5.0.1 Momentum distribution 

In Fig. 1^ the squared 5- wave Green function multiplied by p^ and the 5'-P-interference 
term are displayed for four different energies. The former characterizes the shape of 
the momentum distribution da /dp, the latter is responsible for the forward-backward 
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E [GeV] 


PBorn [GeV] 


iCc [GeV] 


[GeV] 


-3 


15.3 


14.5 


14.5 





20.6 


19.2 


19.3 


2 


26.4 


25.1 


25.3 


5 


34.2 


33.1 


33.3 



Table 1: Approximate position of the peak in the momentum distribution, p^^^^^ is the 
value obtained including rescattering with the full potential, correspondingly pf^sc 
with Coulomb potential, but = 0.187. 



asymmetry and angular dependent spin effects, neglecting of course rescattering. The 
shift towards larger momenta and the narrowing of the distribution with increasing 
energy are clearly visible, as is the strong growth of the P-wave contribution. 

In Fig. ^we display real and imaginary parts of the function if, which is essentially 
given by the ration between S'-P-interference and the squared 5'-wave. In the limit 
of non-interacting quarks this corresponds to p/m^. The function (f^ will become 
relevant in the discussion of angular distributions, the polarization is affected by ip^ 
and 

Fig. ^ shows the influence of the rescattering correction on the momentum distri- 
bution as given in (|73D . The size of rescattering corrections is governed by the three 
functions tpi, tp2 and ip^ which are displayed in Fig. ^. From the deflnition and the 
way they enter the different expressions the following interpretation is possible: ipi 
governs the size of the additional S'-wave amplitude, 1^2 and ip^, respectively, give the 
P- and D-wave amplitudes generated by rescattering^. 

Fig. ^ clearly demonstrates that the attractive force of the tt spectator on the 
decay products from the primary decay leads to a shift towards smaller momenta 
(the peak positions are listed in Table |l|). The dashed line corresponds to the Born 
result. Fig. ^(a), the other two include rescattering, implemented however in different 
ways. The solid line was calculated using the full potential, i.e. the same that appears 
in the determination of the Green functions. For the dotted line, which is hardly 
distinguishable from the solid curve, a pure Coulomb potential with flxed coupling 
constant was used, similar to the previous analyses of [^, |T6[. The calculation for 
the second alternative is somewhat easier because G{p,E) depends on p = |p| only 
and, for the simple choice of the potential, the functions ipi can be converted into one- 
dimensional integrals. Only one integration will have to be performed numerically in 
this case ITB . 



However, to obtain such good an agreement between the dotted and the solid 
curve, a value as = 0.187 has been adopted a posteriori. This large value can easily 
be understood: the scale governing the rescattering corrections should be of the order 



^ Note that their appearance has an origin different from the S"- P-wave interference: the latter 
would vanish for at = whereas rescattering would remain. 
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Figure 6: (a) \pG{p,E)\'^ and (b) the 5'-P-interference term Ipp/m^ • 
Re{G{p,E)F*{p,E)) for rrit = 180 GeV, = 0.125 and three different energies, 
where the sohd hne corresponds to the energy of the IS-Peeik. 
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p[GeV] p[GeV] 

Figure 7: Real (solid) and imaginary (dashed) part of the function ip{p, E) for rrit = 
180 GeV, as — 0.125 and four different energies. The dotted line shows the free 
particle case Re(f — p/rrit. 
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Figure 8: Modification of tlie momentum distribution tlirougli rescattering. Daslied 
line: no rescattering corrections included; Solid line: rescattering contribution with 
full potential included; dotted line: rescattering contribution with pure Coulomb 
potential and ctg = 0.187 included. 
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Figure 9: The four functions V'l, ■j/'r = Re "^i = I^i -^2 and V'3 required for the 
inclusion of the rescattering corrections. 



of top momentum instead of top mass and thus a larger coupling has to be expected. 
The satisfactory agreement between the two schemes is also evident from Figs. and 
[Tl| for the two functions and Note, however, that has been fixed for this 
comparison a posteriori to the value 0.187 which differs from aJ^o?^{M'^m}^) ^ 0.15 
used in [112 



5.0.2 The forward-backward asymmetry 

Ignoring rescattering corrections for the moment the forward-backward asymmetry 
is governed by the relative magnitude of P- versus S-wave amplitudes. The ratio 
between the corresponding two distributions is shown in Fig. |^. For noninteracting 
stable quarks V9(p, E) =\i/mt, and this behaviour persists even in the presence of the 
QCD potential. The energy dependence of (p is not very strong even in the presence 
of the potential, and most of the energy variation of $ visible in Fig. |T2| is driven by 
the energy dependence of the Green function which through the averaging procedure 
samples different momentum regions. 
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Figure 10: The function ipi for two energies using different potentials. Solid line: full 
potential used; dashed line : pure Coulomb potential with fixed as = 0.187 used. 



It should be stressed that the rate and thus the experimental sensitivity varies 
strongly as a function of p (see Fig. |^) and thus only a limited momentum range 
will be explored for any given energy. In fact, given limited statistics only integrated 
quantities will be measured in a first round. These are characterized by the function 



whose real and imaginary parts are shown in Fig. 12. The result depends 
slightly on the cutoff procedure — in a real measurement the cutoff prescription will 
be dictated by experimental considerations. For definiteness all subsequent results 
will be given for the "relativistic" cutoff prescription^ where p/rrit is replaced hj p/E. 

The energy dependence of both the real and the imaginary parts of $ can be 
understood as follows: Above threshold the interaction is of minor importance and 



increases just like p/rrit = \J E/m. It reaches its minimum roughly in the region 
of the would-be- IS-resonance where the 5* wave is still fairly prominent and the 
P-wave is nearly absent. For lower energies $ starts to increase, again roughly oc 



E\/mt, a consequence of the uncertainty principle (see also ||2^). The imaginary 
part of $ which will induce the normal component of the top quark polarization 
levels off at a constant value of around 0.09, well consistent with the perturbative 



prediction (|50| ) $i = 2as(p)/3 with p ^ \J E"^ + Vl. It vanishes rapidly below E = 
real intermediate tt-states are no longer accessible, rescattering is absent and the 
imaginary part vanishes. Throughout most of the threshold region real and imaginary 
part of $ are of comparable magnitude, and so are normal and transverse polarization. 

The ctg dependence of $ is displayed in Fig. |l^. Above threshold the real part is 
practically independent of a^, and entirely determined by the kinematic relation. The 
minimum of the curve is lowered and shifted towards lower values of E, reflecting the 

^We note that this treatment incorporates the bulk of the "relativistic" corrections discussed in 
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Figure 11: The functions and ijj-^ for two energies using different potentials. Solid 
line: ^j^, full potential used; dashed line: ipR^ pure Coulomb potential with fixed 
as = 0.187 used. Dash-dotted line: ipi^ full potential used; dotted line: ipiy pure 
Coulomb potential with fixed ag = 0.187 used. 



decrease of the IS' peak and the larger separation between the "would be" 5*- and 
P-wave resonance. The imaginary part of $ is essentially proportional to as, a result 
obvious from (^Ol). 

The width of the top quark is fixed in the Standard Model: Tf = 1.55 GeV for a 
top mass of 180 GeV New decay modes could increase the width, large CKM 
angles for the mixing between the third and a fourth family could lower the width. 
The sensitivity of $ towards a reduction of Tt is demonstrated in Fig. O. For an 



artificially reduced width the real part of $ is significantly reduced in the resonance 
region. It is evident that both E and contribute to the effective momentum. The 
oscillatory pattern reflects the reappearance of resonance structures. The changes for 
the imaginary part can also be easily interpreted: The suppression of rescattering 
below threshold becomes more pronounced for the reduced width, leading to a reduc- 
tion of $1 for E below —2 GeV. Above this energy rescattering is accessible and even 
enhanced for the case of reduced width: the effective momentum is smaller and the 
QCD coupling correspondingly enhanced. 

Fig. ITBI displays the influence of Higgs exchange on these predictions. The Yukawa 



potential induced by Higgs exchange plus the hard vertex correction according to ||23 



is included in the prediction for the cross section — the prediction for $ contains the 
Yukawa potential only. The modification of the cross section is quite sizeable, the 
influence on angular distribution and top polarization is negligible. 

With these ingredients it is straightforward to predict the forward-backward asym- 
metry ^FB- It depends on the beam polarization x, the cms energy E and the top 
quark momentum p. The predictions without (dashed line) and with (solid line) 
rescattering are displayed in Fig. |16|. The relative size and sign of the rescattering 



29 



0.25 
0.225 
0.2 
0.175 
0.15 
0.125 
0.1 
0.075 
0.05 
0.025 


0.1 



(a) 4>r(E) 



cut-off dependence 



. 1 1 1 1 1 1 1 1 1 1 1 1 1 






X^^^^ 2 m, J 








3 m/4 




m/3 : 







-4 -3 -2-10 1 
(b) 4>,(E) 








-4 -3 -2-10 1 



2 3 4 
E[GeV] 



Figure 12: The integrated function $(£'), (a) real part, (b) imaginary part. The 
different sohd hnes in (a) correspond to cutoffs of = 2mt, rrit, 3mt/4 and m^/S, with 
the highest curve for = 2mt. For the dashed hnes in (a) and (b) the relativistic 
prescription was adopted. 



30 



0.25 



(a) <I>r(E) 








m,= 180 GeV 



-6 -5-4-3-2-101234 
(b) ^,(E) 








a3 = 0.115 
= 0.125 
= 0.135 



-6 -5 -4 -3 -2 



■10 12 3 4 
E [GeV] 



Figure 13: Real and imaginary part of ^{E) for three different values of a. 
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Figure 14: Real and imaginary part of ^{E) for artificially lowered width. t\^^ is the 



standard model width given in [23 . 
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Figure 15: (a) Total cross section and (b) $r = Re$ and $i = Im$ for two different 
values of the Higgsmass Mh and without Higgs (sohd line) . 
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contribution also depend on momentum, energy and beam polarization. Note that 
^FB is affected by both ipi and ip2 — in contrast to the momentum distribution per 
se, which was influenced solely by This is consistent with the interpretation of tpi 
and ip2 as S- and P-wave decay amplitudes mentioned above. 

In the region where most of the rate is concentrated (see Table |I|) rescattering 
corrections can practically be ignored and the qualitative discussion based on ip alone 
is applicable. For precision studies, however, rescattering may become important. 



5.0.3 Top quark polarization 

Three remarks are important regarding the effect of rescattering on the top quark 
polarization. First, what will be discussed in the following is the quantity that can 
be measured by the moments of the lepton spectrum, which should be understood 
as a practical way of defining top quark polarization. Second, an extremely useful 
result is that the normal component Pn is not affected by rescattering and can thus be 
obtained immediately from Figs. § and |^. As a rule of thumb one predicts V^{p, E) ~ 
-0.07sint9. 

Third, the trivial angular dependence of the various spin components arising from 
S'-P-interference seems to be spoilt by rescattering due the ^/'s terms. In the present 
situation, however, t/'s is always multiplied by a numerically quite small factor, i.e. 
for mt = ISOGeV 

mf l + 2y + Sy' 

Since all ipi are of comparable magnitude, the ip^ contribution can be neglected. Hence 
the angular dependence of the longitudinal and transverse polarization essentially 
remains as before. 



In contrast, the influence of V'2 terms on the longitudinal (Fig. O) and transverse 



(Fig. |I8D components of the "polarization vector" is dramatic. The subleading piece 
of V\\, i.e. the coefficient of cos'd, is dominated by rescattering, essentially because 
the S'-P-interference contribution is proportional to the small quantity C|| and even 
vanishes for |x| = 1, a consequence of the normalization condition for the polarization 
vector. This condition does no longer hold once rescattering corrections are taken into 
account by calculating moments of the lepton spectrum (it would of course remain 
true in the spin projection formalism). The dominant piece of V\\ which is given by 
C|i' remains however unaffected and the subleading term which is shown in Fig. |T7| 
vanishes after integration over the top direction. 

The transverse component is also dominated by the rescattering term as it is 
evident from Figs. ^ and |^ that S-P-interference alone would lead to a negative 
result for non-negative x- 

To summarize: Our analysis has shown that the forward-backward asymmetry 
and the normal component of the polarization vector are (relatively) stable against 
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Figure 16: Rcscattcring corrections to the forward-backward asymmetry. Dashed 
hne: S-P-wave interference contribution; sohd hne: full result. The little arrows 
indicate the position of the peak in the momentum distribution. 
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Figure 17: Subleading part of the longitudinal component of the polarization vec- 
tor Vw for 1? = 0. The solid line shows the complete result for unpolarized beams, 
the dotted hne the 5'-P-interference contribution. The dash-dotted line shows the 
complete result for fully polarized beams, where the S-P-wavc contribution vanishes. 
The arrows indicate the position of the peak in the momentum distribution. 
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Figure 18: Transverse component of the polarization vector V±, for d = 7r/2. The 
sohd and dotted hnes correspond to unpolarized beams, the dash-dotted and dotted 
hnes to fully polarized beams. The upper two curves show the full result, the lower 
two the jS-P-interference contribution. The arrows indicate the position of the peak 
in the momentum distribution. 
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rescattering corrections. They determine the top couphngs 03 and 04 via CpB (oT 
possibly some hints on non-standard CP-violating interactions) and both the real 
and imaginary part of the function $ that contains important information on Ft (and 
thus on new physics or the existence of a fourth generation) and the strong coupling 
constant Og. The longitudinal component of the polarization vector is split into two 
parts. The leading piece C^^ is independent of the production angle and measures 
the couplings ai and 02- The subleading angular dependent part and the transverse 
component are small. They are sensitive to rescattering but vanish after angular 
integration. 



A Rescattering corrections 

In this appendix we will give a more detailed description of rescattering using the 
formalism presented in Sect. ^. For simplicity we will restrict ourselves to diagram 
l^a and mention only briefly where differences occur with respect to diagram §b. 

The quantity of interest is the interference contribution between the leading order 
amplitude (with p = |p|, which should not be confused with the four- vector p, and 
correspondingly for k) 



■Gip,E){Gtip,E) + Gtip,E)) (101) 
and the rescattering amplitude 



(27 

Yv{b) ■ Gt{p, E)G{k, E) {Gt{k, E) + Gi{k, E))D^,{k - p)glGF (102) 
to the hadronic tensor 

Hij = JdFS2it;bW+) j dF^2{i]bW')M^M] + c.c. (103) 

Gt and Gi are the nonrelativistic propagators 



£;/2±po + ^rt/2-pV2m4 



and S^ik — p — b) is the b propagator. Neglecting terms of the order [ko — pof and 
(k — p)^ it can be written as 

oY\k — p — b) ^ - 



2Ei{ko -po- n-t ■ (k - p) - te/2E-,) 
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The ko integration may then be performed immediately: closing the contour in the 
lower half-plane and using Coulomb gauge to avoid spurious poles in the gluon prop- 
agator, one only picks up the pole of Gt{k, E) (putting antitop on its mass-shell): 

ko = ^Hrt-E + ky2mt). 

The appearance of K^-^'^ K^D fj,y{k — p) = A^Doo{k —p)gi,Q implies that only Coulomb 
gluon exchange contributes and the replacement glCpD^Qik — p) —iV{\^ — p) is at 
hand. We thus find 

^ (Wn^Y f -., /+(l-7^)A+7.A-/4l-7^)(^+^- m° 
VygJ 7 (2vr)3"^'^ 2E,{-po - ng • (k - p) - zT,/2) 

v{b)Gt{p,E)G{k,E)V{k-p) 

and insert this into the expression for the hadronic tensor. The [bW^) phase space 
integration that arises is trivial due to (p5D, and we arrive at (only displaying one of 
the two interference terms) 



,2 



(27r)3 4 

G*(p, E)G{k, E)V(k-p)G,(p, E) (Gl{p, E)+G*,{p, E) 

r, ^ Tr[ fil - 7^)A+7.A- i>- W ^^-(1 - 7')A-7,A+i^ 



2mtEi ^ + ng ■ (k - p) + ^1^/2 

Now the po integration has become simple: closing the contour in the upper half-plane, 
only the pole of G\{p^ E) is enclosed (putting top on its mass shell): 

po = l-itTt -E + pV2m) Gt{p, E) = ^ . 
2 iTt 

The terms proportional to ko or po in the numerator are thus obviously subleading, 
but the same is also true for the space components. Hence only 

4 7 2Ei{ni ■ {k- p) + tTt) ^ mt 



'■w 



7 ■ (k - p ) 
|k-p| 

remains to be determined, which is relatively straightforward and yields 



(A(k,p)7°- ^,,^_ P^ i?(k,p))(l-y) 



A(k,p) = -i ,,^^ , arctan— — (104) 
|k-p| Tt 

/ |k — p| \ , , 

5 k,p = -K- — -^{l-T. — ^arctan ' ' ) 105 
|k — p| ^ |k — p| 1 1 ^ 
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with y = m^^/mf and /t = (1 — 2y)/{l + 2y). So we finally find the hadronic tensor 
to be 



Hi. 



mtJ (27r)3 
(A(k,p)7° 



G*ip, E)G(k, E)y(k-p)Tr f(l - 7')A+7.A_ 



7 • (k - p) 



1 + 75 /+ 



|k-p| -5(k,p))(l-7')A-7,A+ 2 



plus the corresponding complex conjugate term. The trace can be further simplified 
using / = mf7° which is valid to the order required here, and thus 



A, 



1 + 7V+ / 



m 



(1-7')A+ = A+ 



1 + s+ ■ 77^ 



leading to 



d^k 
(2^ 



xTr 



A+7,(A(k,p 



G*{p,E)G{k,E)V{k-p 
(k - p) ■ 77^ 



5(k,p) 7, 



1 + s+ ■ 7 7^ 



(106) 



|k-p| ^ '^'^ " 2 

Performing the trace, contracting with the leptonic tensors and taking twice the real 
part to include the complex conjugate term, one finally obtains the formula given in 
Sect. ^. The relation between the functions A, B on one hand and ipi, ip2 on the 
other is quite obvious: 



^ , d^k ..^(k,^;) 
2Re /^V(|k-p|)^^-l^A(k,p 



(107) 



(27r)3 

^(|k - Pl)f#^i^5(k, P). (108) 



(27r)3 - |p||k-p| Gip,Ey 

The calculation of the second diagram is essentially the same, the main difference 
arising from the change in the ordering of the vertices which produces another trace: 



Tr 



A_ 



1+75/4 



(A(k,p)7° 



7 ■ (k - p) 



i?(k,p) (l-7')A+7.7, 



2 v----'-^' |k_p| 
and consequently leads to a different contribution to the differential cross section. 



B Calculation of the moments of the lepton spec- 
trum 

It has already been mentioned in Sect. § that for the calculation of the moments of 
the lepton spectrum, the hadronic tensor is replaced by either 

fdpo fdW^ 



Hi 



2ti J 27r 



rfPS2(t;&P^+) JdPS2{W+;lu) JdFS2{t]bW-)^ MiM] (109) 
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or 



(110) 

where the first expression leads to the differential cross section, the second one to the 
moments. The amplitudes M. describe production of the ti pair and its subsequent 
decays: the three-body transition t hlv and the two-body transition i hW~ . 

Let us consider for illustration the simplest case first - the leading vector-vector 
contribution to the hadronic tensors, neglecting terms of order (3. For this contribution 

MM] = glD^{^y\Gip,E)\'\G,ip,E) + G,ip,E)\'T"f' 



X Tr 

where 



7/3 ^a^+l^A- i- 7')A-7,A+(1 + r) (111) 



T"^ = lTr[^7" ^^(1-7')]. (112) 



and the narrow width approximation is employed for the propagator of 

1 



D 



w 



w 



-5{W'-mi^). (113) 



Integration over the two-body phase space PS(t; hW ) has already been described in 
Sect. ^, see (58), and over W"^ is trivial. The po integral can be solved by closing the 



integration contour in the upper half-plane, thus enclosing the two poles 



with the result 



/ -^\Gt{p,E) + G,{v,E)\' = -. (114) 
Using Lorentz symmetry one obtains 

/ dF^2{W] Iv) T"^ = ^{W^W^ - W^g"'^) (115) 

J 24:71 

and 

JdPS2{W-Ju)T'"^ (nl) = ■^[{Wn){W''Wf^ - W^g"'^) - ^Wh^f^P^npWs) (116) 

and integration over the two-body phase space PS2(t; bW~^) also presents no problem 
if one exploits Lorentz symmetry. We have to calculate 

dPS2{t;bW+)[ yV+ l^yV+ + 2m^,y ^) = J dPS2{t;bW+)[ 1/1^1/+ 2mly ^] 
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where 6^ = is assumed, and 



dPS2{t-bW^ 



2 



The calculation of these expressions can be reduced to the calculation of tensor inte- 
grals of the form 



r ■PS(t; hW^) 



(117) 



and 



dPS2(t; hW^)h''W'^ 



12 



2{l + 2y)et^ + {I - y)mig''^ -V^iit^hW^) (118) 



Putting everything together and identifying BR(Vr — > lu) with BR()f: — > hlu), we find 



Hij-n 



2BR(t 

2BR(t blu) 
(4(nt) ^- 



5/.)Hl^Tr 



f)A_7,A+(l + 7') (119) 



|G(p^l^% + %^l 
4(1 + 2|/) 3 



mt 



Tr 



A+7iA_ 



f)A_7.A+(l + 7') 



(120) 



The traces can now be performed. Contractions of the resulting expressions for the 
tensors Hy and Hij.n with the leptonic tensors L*-'^, c.f. (|^), yield the results given 
in Sect. |[ 

Contributions to the moments from rescattering in the tb system, c.f. Fig. |^a, can 
be calculated in a way analogous to that described in App. ^ for the integrals over 
ko, po and PS2(t; bW~) whereas integration over PS3(i:; blu) is performed as described 
in this appendix. It is not surprising therefore that one arrives at final expressions 
which are analogous to those obtained using the formalism of App. ^ The same is 
true for all factorisable contributions, in particular for S-P-wave interference ones. 

A new feature arises for bt rescattering, see Fig. P^, which is non-factorisable. 
Instead of ( |117| , p.l8| ) two new tensor integrals appear in the calculation 



-'I/O 



h la 



dFS2it;bW- 
dFS2{t; bW- 



{b + k — pY 



{b + k — py 



m- 



2 

1-2/ 



t;PS2(t; bW- 
dFS2{t; bW- 



2b{k-p) 



2b{k-p) 



(121) 



(122) 



The fco and Pq integrations can again be performed by closing the integration contours 
in the complex planes. However, in the present case the contours are closed in the 
upper half-plane and the lower p° half-plane, and this choice implies k^—pQ ^ iTt- 
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The integrals /f are proportional to A{k, p) or -B(k, p) again. One finds 

The integrals 12'^ are given by the following formulae: 

7- = 2A(k,p)^il^ (125) 

{k — py{k—p)K mt{l — yY 



|k-p|2 ^ ' "^Tj 647r 
l + 2y X .mt{l-yf 



'-t fr ^fc- - p) - yy x ^1 - -t/J , „^ 

l-2y|kr^t'^-"^ |k-pp i^^^'P%-64;r- ^^2^^ 

The second term in 12 is subleading, because the relevant momentum transfer should 
satisfy |k — p| 3> F^, i.e. the top width only acts as infrared regulator and can be put 
equal to zero where possible (the overall factor l/F^ common to all integrals drops 
out in the final result). The first term must be kept however, and the appearance of 
this second rank tensor is the reason why the function -03 has to be introduced. 
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